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ABSTRACT 

We address two outstanding issues in the sequential accretion scenario for gas 
giant planet formation, the retention of dust grains in the presence of gas drag 
and that of cores despite type I migration. The efficiency of these processes is 
determined by the disk structure. Theoretical models suggest that planets form 
in protostellar disk regions with an inactive neutral "dead zone" near the mid 
plane, sandwiched together by partially ionized surface layers where magnetoro- 
tational instability is active. Due to a transition in the abundance of dust grains, 
the active layer's thickness decreases abruptly near the ice line. Over a range 
of modest accretion rates (~ 10~^ — 1O~^M0 yr~^), the change in the angular 
momentum transfer rate leads to local surface density and pressure distribution 
maxima near the ice line. The azimuthal velocity becomes super-Keplerian and 
the grains accumulate in this transition zone. This barrier locally retains proto- 
planetary cores and enhances the heavy element surface density to the critical 
value needed to initiate efficient gas accretion. It leads to a preferred location and 
epoch of gas giant formation. We simulate and reproduce the observed frequency 
and mass-period distribution of gas giants around solar type stars without hav- 
ing to greatly reduce the type I migration strength. The mass function of the 
short-period planets can be utilized to calibrate the efficiency of type I migration 
and to extrapolate the fraction of stars with habitable terrestrial planets. 
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Subject headings: planetary systems: formation - solar system: formation - stars: 
statics 



Introduction 



In the previous papers of this series Jlda fc Lin|[2004al lbl. boosl . boosl . hereafter Papers 



I, II, III, and IV), we constructed a numerical scheme to simulate the anticipated mass and 
semimajor axis {Mp — a) distribution of planets based on a comprehensive treatment of the 
sequential planet formation scenario. In the scheme, we first generate a set of protoplan- 
etary disk models with various surface densities and depletion timescales, on the basis of 
radio observations. For each disk, we randomly select semimajor axes of the protoplanetary 
seeds and integrate the growth of the protoplanets due to planetesimal accretion, assuming 
that planetesimals have been already formed from dust grains in the disks. If their masses 
become large enough, gas accretion onto the planets is also added. The planets' orbits evolve 
through type I and type II migrations. We integrate growth and orbital evolution of planets 
independently, neglecting dynamical interactions between planets (see discussion in §3.3). 
For the integration, we use the semi-analytical prescriptions based on detailed numerical 
simulations. For details, see §3.3, Papers I and IV. In Paper I, we presented calculations for 
solar-type stars by neglecting the effect of type I migration on the basis that its efficiency is 
poorly determined. With the same assumptions, we simulated the Mp — a distribution for 
stars with a range of metallicity ([Fe/H]) and mass (M^,) in Papers II and III, respectively. 

As a consequence of their tidal interactions with surrounc iing disk gas, embedded em- 



bryos more massive than Mars migrate towards their host stars ( iGoldreich &: Tremaindll980 



Wardlll986l ). In Paper IV, we have considered the influence of this type I migration on the 
planet formation process. Our results indicate that when this effect is fully taken into ac- 
count, the icy cores have a tendency to migrate into their host stars before they acquire 
adequate mass to initiate efficient gas accretion. In order to preserve a sufficient fraction of 
cores which can subsequently evolve into the observed population of gas giants around solar 
type stars, we introduced a "type I migration reduction factor" Ci. The magnitude of Ci 
smaller than unity lengthens the actual magnitude of the type I migration timescale relative 
to that deduced from linear theories (Paper IV; also see the results in §3.3 and 3.4). With 
a range of small Ci ~ 0.03 — 0.1, we were able to simulate a planetary Mp — a distribution 
which is qualitatively consistent with that observed by radial velocity survey. 

While several suppression mechanisms for type I migration under various circumstances 
have been suggested (see references in Paper IV), the extremely small values of Ci we have 
adopted in our previous models remains a challenge to our theoretical construct. In addition. 
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there remain some quantitative discrepancies between the results of our population synthesis 
and the observed a distribution of extrasolar planets. Around solar-type stars, existing data 
show an steep up turn in the frequency of detected planets wit h period around 1-4 yrs 
(corresponding to a in the range of 1-3 AU) (ICumming et al.ll2008l ) whereas the predicted a 
distribution from our past models is essentially logarithmic (see the results in §3.3). 

Another minor issue concerns with the availability of the building-block material for gas 
giants' cores. Current theory of gas giant formation requires a critical mass (Merit) in excess of 
several earth masses in order for the onset of efficient gas accretion (e.g., Paper IV). But, the 
growth of embryos is limited by an isolation mass Mjso. For the minimum mass solar nebula 
(MMSN) model, the surface density distribution of the solid material (S^ oc a~^/^) implies 
Miso is an increasing function of the semimajor axis a. The growth timescale for the embryos 
'i'c.acc also increases with a. On the gas depletion timescale Tdep, the most massive embryos 
emerge near the ice line of a disk comparable to MMSN with a mass Mc ~ Miso < M^it 
(Paper I). The r nagnitude of Msn oc J ^J"^ and the embryos' characteristic growth time scale 



is Tc^acc (IKokubo fc Idall2002l ). In principle, Msn can exceed Mcx \ t in disks which 



are e verywhere several times larger than the MMSN (iKokubo fc Idal 120021 : iThommes et al. 



20071 ). But, the low efficiency of converting these building block material into embryos and 
few giant planets becomes a new issue. The retention of a large amount of heavy elements 
throughout the disk may also lead to noticeable metallicity dispersion araong stars within 
any g iven stellar cluster, contrary to its observed upper limit f Wilden et al. 2002; Shen et aD 
2005h . 



In order to resolve the issues of Mp — a distribution and the adequacy of planet-building 
blocks, we consider an additional effect through which grains and planetary embryos may be 
trapped in some special location in the disk. This effect is associated with an intrinsic local 
pressure maximum in the disk which induces the gas to attain a super-Keplerian azimuthal 
speed. The hydrodynamic drag by the gas on the grains causes them to stall and become 
trapped near the ice line (IKretke &: LinI 120071 ). In §2.1, we describe this grain trapping 
process and show how it may lead to the local (rather than global) enhancement of near 
the ice line such that the formation of gas giants is possible in disks similar to the MMSN. 

In Paper IV, we showed that the formation probability of gas giant planets and hence 
the predicted mass and semimajor axis distributions of extrasolar gas giants are sensitively 
determined by the strength of type I migration. Here we show that the structure of the 
disk near the ice line can also locally stall the type I migration of the embryos. In §2.2, we 
show that the modification of the azimuthal speed of the disk gas near the ice line provides 
another effective mechanism for the retention and rapid accumulation of embryos. In §3, 
we incorporate this effect into our numerical scheme and simulate the mass- semimajor axis 
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distribution of gas giant planets around solar-type stars. In §4, we summarize our results 
and discuss their implications. 



2. Disk structure and heavy elemental retention 



Gas giant planets form in protoplanetary disks prior to the gas depletion. The first 
step in this process is the condensation, growth, and retention of grains. In most regions of 
the disk, a negative pressure gradient in the radial direction reduces the azim uthal speed of 



the g as below the Keplerian value, so the grains undergo orbital decay (e.g., lAdachi et al. 
19761 ). At a ~ 1 AU in a MMSN, the meter-size particles migrate into the Sun within a 



few centuries. Particles can grow faster t han their orbital decay provided the local is 
larger than the surface density of gas Eg (ISupulver fc Linll200ll ). Under these conditions 
and in the absence of turbulence in t he disk, the gr a ins may also settle into a thin layer 
and become gravitationally unstable (jSafronovl 1 19691 : iGoldreich fc WardI Il973l ). However, 
Kelvin - Helmholtz instability caused by velocity d ifference between dust-rich and dust-poor 



layers ( Weidenschilling &: Cuzzj 1993 : Sekiva 1998h an d intrinsic turbulence such as Magneto- 



Rotational Instability (MRI) (IBalbus fc Hawleylll99ll ) may prevent the concentration of dust 
grains. The main challenge in this scenario is how to enhance a local concentration of dust 
grains JVoudin fc Shulbood : Icaraud fc Lin 1 120041 : ICaraudlboOTh . 



If the orbital decay of icy grains due to gas drag is terminated locally near the ice line, 
they would be concentrated near the ice hne at orbital radius aice (eq. [9]). One possible 
mechanism to terminate their decay is through sublimation of icy dust grains and viscous 
diffus ion of the water vapors across the ice line (jStevenson fc Lunind Il988l : ICiesla fc Cuzzi 



20061 ). Across Ojce, an equilibrium may be established in which the outward diffusion of 



water molecules and their recondensation into grains is balanced by the inward orbital decay 
of the solid particles. The exchange of latent heat associated with these phase transition 
significantly reduces the local temperature gradient (Kretke & Lin, in preparation). Here, 
we consider another possible mechanism to terminate the decay, positive radial gradient of 
disk gas due to change in activity of MRI across Oice. This mechanism also halts type I 
migration of planetary embryos and enhances formation and retention rates of gas giants. 



2.1. Layered-accretion and dust retention near the ice line 



The magnitude of gas surface density Tig is determined by the efficiency of angular 
momentum transport in the disk. A leading mechanism for angular momentum transfer in 
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typical astrophysical disks is MRI. However, near the mid plane at a few Alls in protoplan- 
etary disks, the gas may have a sufficiently low ionization fraction which quenches the MR I. 
This region is commonly referred to as the "dead zone" (jGammidll996l : ISano et al.ll2000l ). 
Nevertheless, the surface layers of these disks are exposed to the io nizing photons from the 
central host stars and nearby massive stars (iGlassgold et al.l 119971 ) as well as cosmic rays. 
The ionization fraction x in this layer is sufficient to provide an effective angular momentum 
transfer process which would allow an accretion flow wit h a flux M ~ 10~^-10~ '''M(t^ yr~^. 



which is inferred from observation of T Tauri stars (e.g., iHartmann et al.l 119981 ). But, the 



magnitude of x is determined by an ionization equilibrium which is determined by both 



(Sano et al. 


2000; 


Turner et al. 


2007) 



iterior to Oice, the sublimation of ices 
greatly reduces the effective cross section of the grain population from that outside a^ce- 
Consequently, the active layer is much more extended interior to a^ce than exterior to a^^e- 

The sublimation of ices at r < Oice implies that t he angular momentu m transfer efficiency 
in the surface layers undergo a transition across a\n^ (IKretke &: Linll2007l). In regions of disks 
well interior to the radius of maximum couple (ILynden-Bell &: Pringlelll974l ). the disk flow 
(M) rapidly adjusts to a state of quasi steady state in which M is approximately independent 
of orbital radius r. In principle, the structure of a MRI-driven-turbulent disk requires multi- 
dimensional simulation, even in the limit of steady disk flow. Such simulations are time 



consu ming and are yet to be fully resolved (IFromang fc Papaloizou 1 120071 : iFromang et al 



20071 ). In the present context, the dominant effects on the disk structure can be illustrated 
with the standard ad hoc a prescription to approximate the effective "turbulent viscosity" 
with V = aCsh, wher e c., and h are sound veloc ity and scale height of the disk and a is an 
efficiency parameter (jShakura fc Sunyaevlll973l ). In a steady accretion disk. 



In regions wher e the surface layer is sufficiently ionized to be affected by the MRI turbulence 
(IGammid Il996l ) but the interior is essentially shielded and r e mains inactive (or dead), we 
find it convenient to adopt a prescription by iKretke fc Liru (120071 ) in which the effective 
magnitude of the a parameter. 



a 



SaOa + (S„ - SA)aD 



(2) 



where Sa is the gas surface density of the active layer and and aj) are alpha parameters in 
MRI active and dead zones. Outside the ice line, the condensation of the grains significantly 
reduces the column density of the active layer. Since aA/ao ^ 1, the effective magnitude of 
a decreases as r across the ice line. In a quasi steady state where M is nearly constant of r. 
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Hg (X a ^ so that can be enhanced significantly across the ice fine. This positive gradient 
of can lead to a local pressure maximum. 

When the exponent of disk midplane pressure P, defined by s = cilogP/dlogr, is 
positive, the grain s experience a tail win d and the hydrodynamic drag leads to their outward 
orbital expansion (INakagawa et al.lll986l ). For the equation of state of an ideal gas, P oc pgT, 
where pg and T are mass density and temperature of the disk. Since pg ^ T,g/2h ^ T,gQK/2cs, 
where f^K is Keplerian frequency, positive s is realized when 



p > -q/2 + 3/2, 



(3) 



where p = dlog'Eg/dlogr and q = dlogT / d\ogr (dlogCg/dlogr = q/2). For q = —1/2, this 
condition is p > 7/4. A shallower surface density gradient is required if the release of the 
latent heat essentially suppresses the local temperature gradient. 

Modification of the profile is confined to the ice line region. Since the positive s is 
local, dust grains accumulate near the outer edge of the positive s regio n, which may lead to 
a large enhancement in near the ice line region (IKretke fc Liru 120071 ). Since the isolation 
core mass and the core's growth rate increases with E^, it is possible to build up sufficiently 
massive cores to start the runaway gas accretion process while the residual gas is depleted 
to surface densities comparable to that of the MMSN. 



2.2. Local pressure maximum and type I migration, planetesimal 

accumulation 



Another important consequenc e for the modificati on of the rotation law is the suppres- 
sion of type 1 migration efficiency (IMasset et al.ll2006l ). Even in the stage in which planet 
accretion from planetesimals proceeds, large amount of small dust grains may still fioat in 
the disk ( Tanaka et al. 2005[), so em bryos undergoing type I migration can also be trapped 
near the ice line (jZhang et al.ll2008l ). 



For type I migration, since curvature of the system also affects the locations of Lindblad 
resonances, the condition for outward migration is slightly modulated from eq. The 
migration rate of a planet with mass Mp is given by (ITanaka et al.ll2002l ) 



dv 

— ~ 1.08(p + 0.80g-2.52) 

iJjL 



(4) 



This expression includes contributions from both Lindblad and corotation resonances. The 
dependence on q incorporates uncertainties associated with nonlinear effects. The condition 
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for outward migration is p > — 0.80g + 2.52, which is similar to that for gas drag migration 
(eq. [3]). For q = —1/2, dr/dt > for p > 2.92. As shown below, this condition can be 
satisfied near the ice line. 



3. Simulations including the effect of ice line retention 

3.1. Prescription in the numerical scheme 

The dynamical evolution of the planetesimals is regulated by the surface density distri- 
bution of the gas in the disk. In principle, we should compute the evolution of the gas in 
terms of the standard diffusion equation which takes into account the gas infall, depletion, 
initial and boundary conditions. For computational convenience, we assume steady accretion 
flow and approximate 

M = 3 X 10-V,,oexp(-t/rdep) [M^/yr], (5) 

and calculate evolution of T,g with eq. ([T]), by taking into account the effect of spatial non- 
uniformity of a. T he value of 3 x 10~^Mf7) /yr is typical values of observed M around stars 



with ages ~ Myr (IHartmann et al.lll998l ). In this expression, we assume an exponential 
depletion of the disk gas on some characteristic depletion time scale (rdep). For the purpose 
of illustrating the gross effects of disk evolution on planets' migration, these approximations 
are adequate. However, we have not considered the aspect of non-steady nature of the flow 
which can lead to an outward type II migration at large distances from the host stars. Such 
an effect will be considered in the future. 

Similarly, for computational convenience, we adopt a fiducial prescription for the initial 
surface density distribution of rocky and icy materials as in our previous papers such that 

Srf = 10r/ieeAo(r/lAU)-='/2 1^/^^% (6) 

where the enhancement factor (r/icc) is introduced to take into account the condensation of 
icy grains. In our previous papers , we simply set ?7ice = 1 for r < aice and r/ice = 4.2 for 



r > a-ice according to lHayashil (119811 ). The phase transition between condensed and vaporized 
grains occurs over a small radial range, probably comparable to the vertical density scale 
height h. Here we smooth out the change in r)ice in terms of a tanh function with a width 



h (see Figs. [T]). The enhancement factor may be slightly smaller (~ 3.0) (jPollack et al. 



19941 )]. We found that the retention efficiency of cores against type I migration is reduced 
only slightly and overall features of our results do not change, even if we use ?7ice = 3 for 
r > Oice in our simulation, because a barrier for migration is a more important effect of the 
ice line than the enhancement of S^. 
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In this paper, we consider the potential effect of Eg variation (due to changes in rjicc) 
across the ice hne. In principle, we should take into account any possible accumulation of the 
trapped grains near the ice line prior to this seemingly arbitrary initial state. However, the 
local pressure maximum in the distribution appears only after the gas accretion rate in 
the disk has declined below a critical value. Equation corresponds to the surface density 
of the grains when this barrier at the ice line first appears. Thereafter, as grains congregate 
near the ice line, the ratio of S^/Sg will increase. Eventually, the grain s exert a significant 
torque on the gas to slow down its radial velocity (INakagawa et al.lll986l . Kretke et al. 2008). 
This process further modifies the enhancement near the ice line. 

In a turbulent disk, small (sub-mm) size grains are suspended in disks with scale height 
comparable to h. As approaches to E^, their coUisiona l growth time scale reduces be- 
low their orbital evolution time scale (jSupulver fc Lid 1200 ll ). The vertical velocity shear in 
the disk is also reduced. Then modest (>mm) size grains can sediment from the turbu- 
lent surface layers to the more quiescent midplane regions of the disk and largest particles 
may settle to a suff iciently thin disk to eventually become gravitationally unstable, and 
form planetesimals (Youdin fc Shul l2002l ). against the Kelvin-Helmholtz instability barrier 
( Weidenschilling fc Cuzzi Il993 ). Possible molten surfaces of grains in this region may lead 
to sticky dust collisions, which also enhances dust sedimentation. The critical conditions 
for the onset of gravitational sedimentation and instability of the dust layer depend on the 
magnitude of S^. In those models which we take into account the effect of accumulation of 
icy grains and water vapor near the ice line, we further increase rjice in (Fig. [T]d) by a 
factor of 

/ice = [1 + 2exp(-(r - a,,,)yh% (7) 

As shown later, the enhancement does not affect the results significantly, by the same reason 
of the case of T^ice = 3. We assume that dust grains have surface density that is comparable 
to that of planetesimals, in order to highlight the effect of the ice line. 

Following the simple prescription of our previous papers, we adopt the equihbrium 
temperature in optically thin disk regions (jHayashilll98ll ). 

r \ 1/4 



T = 280 



lAU 



K, 



where L^, and Lq are the stellar and solar luminosity. The ice line is determined by this 
temperature distribution as 

) AU. (9) 



2.7 



Note that the magnitude of a ,ce may be modified by the local viscous dissipat ion (ILecar et al. 



20061 ) and stellar irradiation (IChiang fc Goldreichlll997l : iGaraud fc Linll2007l ). The predicted 
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semimajor axis distribution of extrasolar planets reflects the magnitude of aice- Furthermore, 
Oice changes with time due to disk evolutio n. The movement of the ice line may affect fina l 
mass and orbital configuration of planets (IKennedy et al.ll2006l : Kennedy &: Kenyorul2008l ). 
Inclusion of this potentially important effect is left to a future paper. In the present paper, 
our purpose is to highlight the possibility of halting type I migration near the ice line and a 
great quantitative accuracy is not important. 



3.2. Surface density distribution of the disk gas 

From eqs. ([2]) and ([S]), is given explicitly by 



M 



J] = 

^ ao I lO-^Mo/yr / VSOOK 



lAU 



-g-3/2 



(10) 



where Ti is disk midplane temperature T at lAU. In Figures [T], the equilibrium distri- 
butions are given as a function of M. In this example, we adopt a set of assumed values 



OA 



10 ^, aj) = 10 and Sa where 



mm 



lAU 



[g/cm 



-2i 



This prescription fo r Sa (Sa < Eg) is an order of m agnitude higher than that adopted by 
( iKretke fc Lidl2007l ). Although iKretke fc LinI (120071 ) assumed /xm-size grains, the ionization 
degree in the disk sensit ively depends on grain sizes and t he grain growth significantly 
enlarges the active layer (jSano et al.ll2000l : iTurner et al.l 120071 ). Dispersion in Sa due to the 
range of disk mass is generally smaller than variations in the dust growth properties. For 
the purpose of exploring and highlighting importance of the effect of the ice line, we use this 
high values of Sa. 

Figures [1^ and b show evolution of Hg without and with the enhancement of dust 
grains near the ice line. In the regions where ^ Sa, the disk is mostly MRI dead and 
a ~ while a ~ oa in the S^ ~ Sa regions. During early stages of disk evolution (epochs 
with high M), the active surface layer occupies an insignificant fraction of the total disk, 
so the transition across aice does not modify the disk structure nor produce positive radial 
gradient. At advanced stages of disk evolution when Sg is very low (M < 10~^Mq yr~^), 
the active layer extends throughout the entire extent normal to the plane of the disk and the 
magnitude of S^ is also a monotonically decreasing function of r. However, for intermediate 
stages (1O~^M0 yr^^ ^ M < lO^^M© yr"^), the S^ distribution has a local maximum with 
positive gradient near the ice line. In the case with the S^ (/ice) enhancement, the bump in 
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Tig is more pronounced. The power index for the Tig gradient is more than 2.9 near the ice 
hne. Based on the consideration presented in §2.2, it is suggestive that type I migration of 
the planetesimals may be halted near there. 

Note that the trapping efficiency of both grains and cores is determined only by the value 
of p and it does not depend on the amplitude of the surface density variation near the ice line. 
The most effective trapping location is interior to the ice line where the pressure gradient is 
most positive. The actual range of M for effective ice-line barrier is a function of the poorly 
understood od- Furthermore, the size distributions of dust grains due to coagulation and 
fragmentation also regulate the magnitude of T\ and when the positive pressure gradient 
appears, but the evolution of the size distributions is also poorly understood. Nevertheless, 
the qualitative implication of this physical effect is unavoidable, the local pressure maximum 
near the ice line is essential for the retention of the planet-building blocks. 

Figures [2] show the time evolution of planetesimals' mass and semimajor axis for the 
model with fg^ = 3 (without the enhancement). For illustration purposes, the planetes- 
imals' initial semimajor axis distribution is chosen to be separated with equal logarithmic 
intervals. This figure shows that some cores are trapped near the ice line and grow there 
until they start to undergo efficient gas accretion. 

3.3. Population synthesis 

With these prescriptions, we carried out a series of Monte Carlo simulations to generate 
theoretical mass and semimajor axis distribution for extrasolar planets. Except for the disk 
gas surface density Tg and the q dependence of type I migration rate (see below), the formulae 
for growth and migration are identical to those in Paper IV. 

We first generate a 1,000 set of disk models with various values of fg^Q (the initial value 
of fg) and Tdcp. In this paper, we focus on solar-type stars and consider a log uniform 
distribution for in the range of 0.8-1.25Mq. We assume = Lq{M^:/ MqY and fg^ oc 
Ml (Paper III). These choices of M^, dependences do not significantly affect the results, 
because the range of for these models is relatively narrow. The luminosity also 
changes with time during pre-main sequence phase, in which planet formation proceeds. It 
affects the location of Oice, but a\ce is affected more by disk structure, as mentioned in §3.1. 
So, we neglect evolution of L^, for simplicity. 

We also assume that raop has a log uniform distribution in the range of 1-10 Myrs. 
This assumption is based on the observation that in young clusters within this age range, 
a fraction of coeval stars have signatures of disks and this fraction declines linearly with 



the cluster age (IHaisch et al.ll200ll ). The modest dispersion in this fraction for clusters with 
similar ages is another indication that there is a considerable spread in the magnitude of T^ep 
even in star forming regions with similar initial conditions. 



We adopt the same prescriptions for the distributions of fdfl and fg^ as those in Papers 
II-IV. For the gaseous component, we assume fg^ has a log normal distribution which is 
centered on the value of fg^ = {M^/MqY with a dispersion of 1 ((51og^o/5,o = 1-0) and 
upper cut-off at fg^ = 30, independent of the stellar metallicity. For the heavy elements, 
we choose fdfl = lO^^^^^^^/gfl, where [Fe/Hj^ is metallicity of the disk. We assume these 
disks have the same metallicity as their host stars. In this prescription, throughout the 
entire disk varies with [Fe/H] in contrast to the local enhancement near the ice line in some 
enriched models. 

For each disk, 15 semimajor axes of the protoplanetary seeds are selected from a log 
uniform distribution in the ranges of 0.05-50AU. The averaged orbital separation between 
planets is 0.2 in log scale, or equivalently, the averaged ratio of semimajor axes of adjacent 
seed planets is ^ 1.6. 

For each set of fgfl, fdfl and a values, we integrate the protoplanets' growth through 
planetesimal and gas accretions with semi-analytical prescriptions based on detailed numer- 
ical simulations. The initial mass is arbitrarily set to be a small value, Mc = 10^°g, albeit 
the orbital evolution of planets with such initial mass would not be affected by the effect 
of hydrodynamic drag. The choice of this initial mass does not affect the results because 
the accretion time scale for embryos with masses > 10^°g increases with Mc, i.e., they have 
entered the oligarchic growth phase. 

In outer regions (> lOAU), core growth is slow and scattering of planetesimals by the 
cores rather results in ejection from the systems, so in situ formation of giant planets is 
limited (Paper I). On the other hand, in inner regions, limited amount of planetesimals in 
feeding zones suppresses formation of giant planets (Paper I). Here, individual planets are 
integrated independently, neglecting dynamical interactions between them. Although the 
change in semimajor axis is not significant for inward scattering due to energy conservation, 
it can be very large for outward scattering. So, giant planets in outer regions would have 
been scattered from inner regions and they would have relatively large orbital eccentricities 
unless adequate damping has existed. Since the scattering process is not included in the 
current simulations, predicted distribution in outer regions would have uncertainty in outer 
regions. This should be cautioned when our current results are used for consideration of 
astrometric or direct imaging observations that are sensitive to planets in outer regions. 
The effect of the uncertainty is smaller for radial velocity observations that are sensitive 
to short-period planets, because these short-period planets would have undergone much 
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greater orbital migrations due to tidal interactions with disk gas (that are included in our 
simulations) than those due to the dynamical scattering. We will include the effects of 
dynamical interactions in a subsequent paper. This uncertainty does not affect the purpose 
of the present paper, highlighting the possibility of halting type I migration near the ice line. 

We integrate the planets' orbital evolution through type I and type II migrations in a 
disk with r^cp around a star with M^. The magnitude of fd at a given location r continuously 
decreases with time from its initial value /d,o as planetesimals are accreted by embryos which 
in term undergo orbital decay. Note that here the semimajor axis a is identified as orbital 
radius r, because we neglect evolution of orbital eccentricities. For the gas component, we 
adopt a prescription for an exponential decay with decay constant Tdep as eq. ([5]). 

While the type II migration speed is determined by an analytical formula with an 
empirical numerical factor (Paper IV), the type I migration speed is given with a scaling 
factor Ci by 

r = Ciriinear, (12) 

where runear is given by eq. (jlj) with q = —1/2. The parameter p in the equation is calculated 
by Eg (eq. [lU]) at each radius and each time. If the calculated p is larger than 2.92, the 
migration is outward. We regards Ci as a parameter and do simulations for different values 
of Ci. We artificially terminate type I and II migrations near disk inner edge where the 
orbital period of the planetesimals is 2 day (~ 0.03AU for M^, = IMq) in a similar manner 
as in Papers I- IV. (The fate of close-in planets will be examined in future investigations.) 

In the present series of simulations, if a planet arrives at the inner edge, a next-generation 
planetary seed would be introduced. The evolution of this new planet would be integrated 
with the residual planetesimal surface density (see Paper IV) . For many models, it is possible 
for several planets to reach the proximity of their host stars prior to the severe depletion 
of the disk gas. In all models, we record not only the individual close-in planets but also 
consider the limiting possibility that all the short-period planets around common host stars 
may undergo dynamical instability after the gas depletion, collide, and coagulate into a single 
entity (see §3.7). 

3.4. Mass and semimajor axis distributions 

In equation f[T^ . Ci represents an reduction factor. For exploration purpose, we consider 
models with Ci in the range of 0-1. The predicted Mp — a distributions of extrasolar planets 
is shown in Figures [3] for Ci = 0.03,0.1,0.3, and 1. Figures [3^ are the results of models in 
which we have neglected the ice-line barrier. In the absence of a bump near the ice line. 
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we assume = 750/g,o exp(— t/rdcp)(r/lAU) ^[g/cm^] and a = 10 ^ throughout the disk. 

In contrast, the simulated resuhs for models that explicitly include the ice line barrier 
are shown in Fig. [Sb (without the enhancement) and c (with the enhancement by a 
factor of /ice as eq. [7]). We adopt [Fe/H]= 0.1 for all these models because the on-going 
radial velocity surveys have been focusing on relatively metal-rich stars. ([Fe/H]-dependence 
is shown in §3.6.) To compare with the theoretical prediction and the data of radial velocity 
observations, we also plot the data of all planets (discovered by radial velocity surveys), 
around stars with M,, = 0.8-1.25Mq. In order to correct for the projection effect, we 
adopted 4/7r (~ 1.27) for the 1/ sinz factor. 

First, we discuss the simulated results in which the ice line barrier has been neglected 
(Fig. [3^1^). These results are essentially the same as those in Paper IV, except that in our 
previous simulations, we adopted = IMq without any dispersion. In the absence of any 
ice line barrier, only for Ci < 0.03 {i.e. in the inefficient type I migration limit), the predicted 
population of gas giants matches well with the observed data (see below). In models with 
higher values of Ci(> 0.3), only the low-mass (Mp < lOM®) cores can survive the type 
I migration. These low-mass cores cannot evolve into gas giants because their envelope 
contraction time scales are generally much longer than the gas depletion time scales (Paper 
IV). 

We now consider the influence of the ice-line barrier. In a series of models with a S^, 
bump (Figs. [3)d and c), much larger populations of gas giants are generated than the barrier- 
free models with corresponding values of Ci (Fig. [3K)- Furthermore, in these models, the 
semimajor axis distribution of gas giants shows a ramp up in the planetary population with 
a ~ 1-3 AU (Figs. [3)d and c). Although in the limit of small Ci, gas giants can also form 
without the ice line barrier, their a distribution is essentially logarithmic (Figs. [3^, also see 
Paper I). 

This dichotomy can be attributed to the role of the ice-line barrier which is to preserve 
cores until they become sufficiently massive (with Mc ~ lOM®) to efficiently accrete gas. 
Planetesimals formed at large distances from their host stars migrate inward and become 
trapped near the ice line. They accumulate and coagulate into embryos at this barrier 
until either the disk gas is depleted or they attain a critical mass to initiate efficient gas 
accretion. However, planetesimals formed interior to the ice line cannot halt their inward 
type I migration until they have reached the proximity of their host stars. If these close-in 
planets are embedded in residual gas, their eccentricities would be effectively damped and 
their further growth would be inhibited by dynamical isolation. After the severe depletion 
of the disk gas, dynamical instability can lead to orbit crossing, cohesive collisions, and 
emergence of relatively massive short-period planets. Although some of these cores may have 
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Mc > IOM0, the depletion of gas would prevent the build up of their envelopes. Around 
solar type stars, it is difficult for gas giants to form in close-in orbits. 

After they have acquired planet masses comparable to those of Jupiter and Saturn, gas 
giants formed near the ice line barrier open gaps and undergo type II migration (Figures [2]). 
However, since the ice line barrier is only effective for a range of modest M ~ 10~^ — 1Q~^Mq 
yr~^, most gas giants emerge in relatively passive disks and the extent of their migration 
is somewhat limited. Consequently, the fraction of stars with short-period gas giants {rjsj) 
is much reduced in the models with the ice-line barrier than the low-Ci models in which 
the ice-line barrier effect is neglected (see further discussions in §3.7). Since the magnitude 
of Vsj predicted with these new models is comparable to its observed value (~ 1%), we no 
longer need to invoke extensive disruption mechanisms for the excess predicted short-period 
gas giants. The reduction in the gas giants' type II migration, due to their late emergence, 
also helps to preserve the upturn in their a (or P) distribution which is a signature of their 
formative environment rather than some evolutionary outcomes. Overall, these simulated 
results are in a better agreement with the observed data. 

There are also many cores trapped near the ice-line barrier which do not have sufficient 
mass to accrete gas prior to its severe depletion. When M reduces below ~ lO~^M0yr~^, 
these cores resume their type I migration due to the removal of the ice line barrier. However, 
the magnitude of of the residual gas is generally too low for the initially trapped cores 
with Mc ~ a few to undergo extensive migration. This modest population of rocky 
planets reside close to the habitable zone (hereafter HZ) slightly interior (at 1-3 AU) to 
the ice line barrier region (see further discussions in §3.7). Since these cores form in nearly 
isothermal regions where gas-phase- vapors and solid-phase-ice coexist, they may contain a 
substantial amount of life-supporting water. The orbits of a fraction of these rocky planets 
may be destabilized by coexisting gas giants in the same systems. But around host stars 
without any gas giants, their orbital eccentricity may be damped by their interaction with 
the residual planetesimals. 

Despite the production of occasional intermediate-mass (Mc ~ a few M^) planets at < 
a few AU, a characteristic "planet desert" (see Paper I) is still prominent in Figures (Sb and 
c. The boundaries of the desert domain are modified by the efficiency of the barrier. Since 
the ice- line barrier cannot act to preserve them, low-mass (Mc < M@) planetesimals and 
embryos formed interior to it migrate to the proximity of their host stars in the presence of 
relatively small amount of residual disk gas. 

We now consider the influence of initial Tj^ enhancement near the ice-line barrier. As 
we stated in the introduction, Miso oc and Xc^acc oc . In principle, the enhancement 
of Ed by a factor of /ice near the ice line leads to rapid growth of embryos with relatively 
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large isolation masses. In the comparisons between these two series of models, we find that 
the distributions of gas giants are similar, albeit a /ice enhancement leads to the emergence 
of marginally more detectable planets, especially in the limit of modest values of Ci(> 0.3). 
The weak dependence is also partially caused by the ice-line barrier being effective only after 
M < 10~^Mq yr~^ so that relatively massive embryos which emerged during earlier epochs of 
active disk evolution would not be preserved. This self-regulated retention condition reduces 
the advantage of the initially more massive disks. The principle factor for the formation 
probability of gas giants is the existence of the ice line barrier rather than any localized 
initial concentration of S^. 



3.5. Enhancement of formation probability of gas giants and domain of 

planetary desert 

In order to quantitatively compare with observations, we determine the fraction {tjj) of 
stars which harbor plari ets within the present radial- velocity detectability limit. Following 



Fischer &: Valentil (120051 ) . we take a conservative estimate on the magnitude of radial velocity 
{vr > lOm/s) and orbital periods {Tk < 4 years) for the detectable conditions. Because we 
artificially terminate type I and II migrations near disk inner edge and we have not specified 
a survival criterion for the close-in planets, we exclude close-in planets with a < 0.05AU in 
the evaluation of rjj. (see further discussions on the short-period planets in §3.7). 

The values of rjj are plotted as a function of Ci in Figure H^. The dependence of r^j 
on Ci is much weaker among models which take into account the ice line effect than those 
without it. Even for C^ = 0.3-1, the predicted rjj can be comparable to the observed data 



( iFischer fc Valentil l2005l ). Migrating cores are trapped near the ice line, almost independent 
of the strength of type I migration, and they accrete planetesimals until either 1) they can 
start runaway gas accretion or 2) Eg near the ice line is reduced below and they resume 
their type I migration. The coupling effect of MRI activity and the ice line barrier can 
enhance formation and retention rates of gas giants without significant reduction of the 
strength of type I migration. 

A close inspection of Figures [2] shows that although the paucity of intermediate-mass 
intermediate-period planets is maintained, the evacuation of this planetary desert is less 
severe in the series of models with the ice- line barrier than those which neglected this effect. 
The upper Mp boundary of the desert shifts downward with increasing Ci and /d,o where a 
large a limit extends to the HZ. A few planets emerge in this sparsely populated region and 
some of them have masses and orbits similar to those of the Earth. In order to extrapolate 
the probability of finding habitable planets from the known extrasolar planets, we adopt the 
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Terre strial Planet Finder's convention for the domains of habitable planets (jPes Marais et al. 



20021 ). We plot in Figure the fraction of stars with habitable planets [r]^) which is defined 
as planets with Mp ~ 0.3 — lOM^ and a ~ 0.75 — 1.8 AU around G dwarfs. 

These planets are formed in the ice line barrier region of those disks which contain a 
marginal amount of heavy elements. Their progenitors' growth toward super-earth-mass 
cores is slow. When they finally acquired sufficient masses to efficiently accrete gas, the 
magnitude of T,g near their orbits also decreases below such that the residual disk does 
not have adequate mass to promote these protoplanets to attain their full growth potential 
prior to the severe gas depletion. However, a relatively small is needed to induce the super- 
earth-mass cores to undergo fractional orbital decay. The late "leakage" of a few cores into 
the planetary desert does not significantly alter the landscape of the planetary desert, albeit 
the presence of a few intermediate-mass intermediate-period planets is no longer forbidden 
(also see previous subsection). 



3.6. Metallicity dependence 

Based on the assumption that protostellar disks attain the same fraction of heavy ele- 
ments as that contained in their host stars, we showed, in Paper II, that gas giant planet 
formation is more prolific around metal-rich than metal-poor stars. In addition to the avail- 
ability of planet building blocks, disks with enhanced metallicities also contain a greater 
population of grains which can modify the ionization fraction, the depth of the active layer, 
and potentially the retention efficiency of the ice-line barrier. But the range of M over which 
the ice line barrier is effective does not depend explicitly on the magnitude of the metallicity. 

Figure [5] shows the values of rjj which are plo tted as a function o f [Fe/ H]. Open circles 



with error bars are observational data compiled by lFischer &: Valentil (120051 ). Triangles with 
dotted line shows the results without type I migration, which are based on similar results 
to those shown by Paper IV. The correlation reflects the fact that higher [Fe/H] enhances 
formation of cores large enough for runaway gas accretion through higher q- In Paper IV, 
we suggested that the effect of type I migration enhances the 77j-[Fe/H] correlation, because 
type I migration is more efficient in metal-poor disks. Squares with dashed line represent 
the results from models in which type I migration with Ci = 0.03 is included and the ice 
line barrier is neglected. A modest rjj ar nplitude and its steep [F e/H] dependence are in a 
better agreement with the observed data ( Fischer &: Valenti 20051 ) than the models in which 



the effects of type I migration has been altogether neglected (represented by triangles and 
dotted lines). 
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We now consider models in which the ice-hne barrier is taken into account. The filled 
pentagons and circles with solid lines respectively represent models without and with the 
Erf o enhancement at the ice line. Despite an order of magnitude increase in the value of 
Ci{= 0.3), the magnitude of rjj for the ice-line barrier models is comparable to or larger 
than those in which the ice line barrier effect is neglected. 

We note from the comparison of the pentagon and fiUed-circle data points that the 
magnitude of r]j increases E^^o by a similar factor for all metallicity. This trend is another 
indication that the ice-line barrier essentially suppresses the depletion by type I migration 
effect once the gas surface density decreases below S^. The greater availability of planet- 
building blocks enhances the probability of gas giant formation in these models. 

3.7. Short-period planets 

In paper IV, we carried out a survey of short-period planets which can migrate to within 
0.03 AU from their host stars. In Figures^ we plot mass functions of short-period planets in 
the cases (a) without the ice- line barrier, (b) with the ice-line barrier and no enhancement, 
and (c) with the ice-line barrier and the E^ enhancement. In the upper panels, we assume 
that all the short-period planets survive without any further coagulation near their host 
stars. After the severe depletion of the disk gas, however, djTiamical instability may develop 
and induce the coagulation of planets. In the lower panels, we consider the limiting case 
that around individual host stars, all planetesimals migrate to their proximity would merge 
into one short-period planet. 

If we neglect the ice-barrier (the upper panel in Figure Efe'), the mass distribution of these 
short-period planets would be biased toward Jupiter mass for Ci < 0.001 (filled circles). 
However, even with a relatively small but finite value of Ci(> 0.01), there would be a 
population of Earth-mass planets. For Ci ~ 0.01, the resultant mass function of the short- 
period planets would be bimodal. For Ci > 0.1, formation of gas giants is inhibited so that 
the mass function is highly peaked at a few M^'s. 

This dichotomy is clearly demonstrated in the upper panel in Figure [7^ where the short- 
period planets are classified into three populations. The filled circles, squares, and triangles 
represent the fraction of stars (77^) which bear short-period planets with Mp in the range 
of > lOOM®, 20 — IOOM0, and 1 — 20M^, respectively. The population of planets with 
Mp > IOOM0 rapidly decreases with Ci, while that of smaller planets decreases more slowly. 
The observed value of rjsj = rjs{Mp > lOOM©) ~ 1% can be matched with Ci ~ 0.1-0.3, 
provided most short-period planets are retained. 
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In these models, Ci < 0.03 is needed to match rjj of all the gas giant planets and the 
over-prediction of 1]^ can still match its observed value under the assumption that more than 
95% of the short-period gas giants have perished. There are also a similar population of 
intermediate mass (10 — lOOM^) and a much larger population of earth-mass short-period 
planets. If all the planets which migrated to the proximity of their individual host stars can 
subsequently collide and merge, the intermediate and low mass planets would become the 
building blocks of the massive short-period planets (see the lower panels in Figures [6^ and 
[7^). But such pla nets are mostly c omposed of rocky cores rather than an extended envelope 
and HD 149026b (ISato et al.ll2005l ) is a good candidate for such massive mostly rocky short- 
period planets. In order to match the observed value of risj, most of these planets must be 
lost to their host stars. 

We now consider the statistical outcome of short-period planet formation in the limit 
that type I migration is suppressed by the ice-line barrier. This process not only promotes 
the emergence of Jupiter-like gas giants but also introduces a strong modification in the mass 
function of the short-period planets (see Figs. [6)d and c). The effect of the enhancement 
is small in the function. Compared with Fig. [6ti, the dependence on Ci is very weak with the 
ice-line barrier. The magnitude of the observed value of rjj can be matched with Ci ~ 0.1 — 1 
(see Fig. and the middle panels of Figs. [3]d and c). With this range of parameter, the 
model values of rjg for gas giants (circles) are also comparable to that observed, provided 
some of them have coagulated and most of them have been retained in the proximity of their 
host stars. The magnitude of r]s{Mp) (Figs. [7)d and c) and the mass distribution of the short- 
period planets (Figs. [Hb and c) indicate that there is a comparable or larger population of 
short-period intermediate- mass (20— lOOM®) planets as that of short-period gas giants. Since 
the value of Mp for the peak of the mass distribution is a function of Ci, its determination 
can provide valuable calibration for the theoretical models. 

The retention efficiency of the short-period planets remains uncertain. Several halt- 
ing mechanisms have been proposed for the short-period planets. They include planet-star 



Gu et al. 


2003; 


Masset et al. 


2006 



Laine fc Linll2008l ). Although the efficiency of these processes is to be determined in the 
future, the lower mass rocky planets are more likely to survive than the gas giants in most 
of these scenarios. Thus, we anticipate the prolific discovery of intermediate to low mass 
close-in planets in the near future searches. We predict the fraction of stars which bear 



short-period earth-mass {rjsE = Vs{Mp 



2OM0)) is comparable to that of all gas giants 



(?7j). This result is insensitive to the initial surface density distribution of the heavy elements 
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Based on these results and the cahbration of rjj for all known planets, we extrapolate 
770 ~ 5-20 % depending on the magnitude of Ci. The magnitude of 77® can be reduced by 
perturbations from gas giants that we have not taken into account in our analysis. However, 
as stated before, for relatively large Ci, the reduction would be small. 

4. Summary and Discussions 

In this paper, we confronted two main outstanding challenges to the core accretion 
scenario: 1) the efficient utilization of planets' building block material and 2) the retention 
of cores to form gas giant planets despite the hydrodynamic drag by the disk gas on grains 
and the orbital decay of the protoplanctary cores due to their tidal interaction with the 
disk gas (type 1 migration). We propose that the resolution of these conundrums lies in the 
existence of an ice line barrier. 

In our scenario, we identify the ice line as a transition zone for the efficiency of angular 
momentum transfer in the gas disk. The presence of additional grains outside the ice line 
reduces the thickness of the active layer and the efficiency of angular momentum transfer. 
Within a range of the disk mass accretion rate (M comparable to that observed for the 
classical T Tauri stars) , these effects lead to local maxima in the pressure and surface density 
distribution in the disk. In the region slightly interior to ice line (closer to the star), a 
super-Keplerian gas flow provides a barrier to retain the heavy elements. Consequently, 
protoplanctary embryos coagulate and grow in mass at this location. Up on acquiring several 
Mq, these cores initiate efficient gas accretion and evolve into gas giants. Since this ice 
line barrier is effective over a range of modest mass transfer rate (10^*^ — IO^^Mq yr~^), 
it regulates the formation of the first-generation gas giants during a preferred epoch at a 
preferred location around their host stars. 

In relatively massive disks, the emergence of the first-generation planet leads to gap 
formation. Just beyond the outer edge of gap, another pressure maximum is induced. Sim- 
ilar to the ice line, this local pressure maximum promotes the congregation, coagulation, 
and rapid growth of later generations of gas giants. However, the physical extent of the 
planetary system is limited to the region where the local Keplerian speed is comparable to 
the surface escape speed of the critical mass (a few M®) cores. Outside this location, close 
(non coUisional) encounters by these cores lead to ejection of their neighboring planetesimals 
from the system so that the depletion of the building block material is faster than the cores' 
growth. 

Although gas giant formation is inhibited in relatively low-mass disks, they may nonethe- 



-20- 



less be the birth places of Earth-mass rocky planets in habitable zones as well as in the 
proximity of their host stars. The final assemble of these planets may proceed throughout 
the age of their host stars. While, in the post T Tauri phase, there is still small amount of 
residual gas in the disk, many of these entities may undergo type 1 migration over limited 
radial range towards their host stars. After the disk gas is severely depleted, these planets 
destabilize each other's orbits, albeit the eccentricities of the most massive planets may be 
suppressed by their gravitational interaction with the less massive residual planetesimals. 

Through physical collisions, a population of a few-Me planets is expected to emerge in 
the habitable zone around 5-20 % of solar type stars. Some close encounters can also lead to 
limited migration. Nevertheless, a desert in the Mp — a distribution is likely to be preserved. 
The paucity of planets with intermediate mass and semi major axis can be established with 
various observational surveys in the next few years. The boundaries of this domain contain 
valuable information on the efficiency of various process and valuable calibrations for the 
theoretical models. 

Many planetesimals and embryos migrate to the proximity of their host stars, especially 
in the limit of relatively efficient type I migration. Many cohesive collisions may also occur 
among these planetesimals and embryos in the proximity of their host stars. Provided most of 

these planets arc retained (despite the effect of their interaction with the magnctosphcre and 
tide of their host stars), we anticipate the fraction of solar type stars with short-period earth- 
mass objects to be several times that of short-period gas giants. The physical composition 
of these close-in planets is most likely to be rocky to icy. Their observed mass distribution 
will enable us to calibrate the magnitude of the type I migration (the magnitude of Ci). The 
next step is to extrapolate and to calibrate the fraction of low-mass M-dwarfs and high-mass 
G sub-giant stars which may bear potentially habitable planets. 
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Fig. 1. — The evolution of equilibrium distributions with the coupling effect of MRI 
and the ice line. The distributions (a) without the effect of enhancement in due to 
the grain trapping and (b) with the effect. The solid lines express the distributions for 
M = 10-^ 3 X 10-^ 10-^ 3 x 10~^ 10-^Mq/jy from top to bottom. The dotted and dashed 
lines are Sa and with = 1. 
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Fig. 2. — The time evolution of planetesimals' mass (Mp) and semimajor axis (a) for the 
case of Fig. [H Here, fg^o = 3 and r^ep = 3 x 10^ years. The bottom panel shows time 
evolution of M. 
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Fig. 3. — The mass and semimajor axis distribution of extrasolar planets. Units of the mass 
(Mp) and semimajor axis (a) are earth mass (Mq) and AU. (a) The results in disks without 
the Eg bump due to the coupling effect of MRI activity and the ice line, (b) those with the 
bump in but without the enhancement, and (c) those with both the effects. The top 
panels are observational data of extrasolar planets (based on data in http:/ /exoplanet.eu/") 
around stars with = O.8-1.25M0 that were detected by the radial velocity surveys. The 
determined MpSini is multiplied by 1/ (sinz) = 4/7r ~ 1.27, assuming random orientation of 
planetary orbital planes. The other panels are theoretical predictions with M^, = O.8-1.25M0 
for various values of Ci. The dashed lines express observational limit with radial- velocity 
measurement precision of Vr = lOm/s. In these models, the magnitude of the metallicity 
[Fe/H] = 0.1. 
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Fig. 4. — The theoretically predicted (a) rjj and (b) rjE (r^e) functions of the type I mi- 
gration strength Ci. rjj is the fraction of stars harboring the planets within the detectability 
limit of by the magnitude of radial velocity {vr > lOm/s) and orbital periods (Tk < 4 years). 
rjE {rj^) is the fraction of stars harboring "habitable" planets with Mp ~ 0.3 — lOM^ and 
a ^ 0.75 — 1.8 AU around G dwarfs. Filled squares with dashed line, pentagons with solid 
line, and circles with solid line express the values for the distributions of Figures [3fe', b and 
c, respectively. (The results for Ci = 0.01 are also included.) 
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Fig. 5. — Metallicity dependence of rjj. Open circles are observational data by 



Fischer &: Valentil (120051 ) . in which bars express statistical errors. The others are theoretical 
predictions. Filled triangles with dotted line and squares with dashed line express the results 
without the Eg bump near the ice line. The former and latter are the cases with Ci = 
and 0.03, respectively. Filled pentagons and circles with solid line express the results with 
the Eg bump. The former and latter are the cases without and with the enhancement, 
in which Ci = 0.3 is used. 
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Fig. 6. — Mass functions of short-period planets in the cases (a) without the ice-hne barrier, 
(b) with the ice-hne barrier and no Hd enhancement, and (c) with the ice-hne barrier and 
the Erf enhancement, which are calculated from the distributions of Figures [3|i, b and c, 
respectively. The filled circles, open circles, squares, and triangles express the results for 
Ci = 0.001,0.01,0.1 and 1. In the upper panels, it is assumed that all the short-period 
planets survive without any further coagulation. In the lower panels, it is assumed that all 
the short-period planets in the same disk coagulate into one body. 
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Fig. 7. — Population of giant (> lOOM^), intermediate-mass (20 — lOOM^), and lower-mass 
(1 — 20M^) planets as functions of Ci. They are expressed by circles, squares and triangles, 
respectively. The meanings of (a), (b) and (c) are the same as those in Figures [61 



